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Most work on wireless network throughput ignores the temporal correlation inherent to wireless channels 



because it degrades tractability. To better model and quantify the temporal variations of wireless network 
throughput, this paper introduces a metric termed ergodic transmission capacity (ETC), which includes spatial 
and temporal ergodicity. All transmitters in the network form a homogeneous Poisson point process and 
all channels are modeled by a finite state Markov chain. The bounds on outage probability and ETC are 



c/2 ' characterized, and their scaling behaviors for a sparse and dense network are discussed. From these results, 

o ; 

we show that the ETC can be characterized by the inner product of the channel-state related vector and the 
invariant probability vector of the Markov chain. This indicates that channel-aware opportunistic transmission 
does not always increase ETC. Finally, we look at outage probability with interference management from a 
' stochastic geometry point of view. The improved bounds on outage probability and ETC due to interference 

management are characterized and they provide some useful insights on how to effectively manage interference 

CO ! 

, in sparse and dense networks. 



I. Introduction 

X 

In the past decade, our understanding of large wireless network capacity has increased considerably, 
but perhaps still comprises more questions than answers, especially for realistic models. Gupta and 
Kumar's landmark work [lj, for example, introduced the transport capacity metric and derived scaling 
laws on it in a size-limited network. Another more recent example is transmission capacity proposed 
in B23 which is a spatial throughput metric for Poisson-distributed transmitters in an infinite network 
with outage constraints. Almost all of the studies following the aforementioned approaches did 
not consider temporal affections. For a wireless network with long-term time-varying channels, its 
snapshot throughput may not provide a full picture of how the throughput evolves over time. 
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In this paper, we introduce a metric capable of characterizing the network throughput induced 
by channels with temporal and spatial ergodicity. This metric is called the ergodic transmission 
capacity (ETC), and it measures the maximum long-term average rate (in bps/Hz) that can be sent 
per unit area in the network with an outage constraint. We evaluate the ETC under assuming all 
transmitters form a homogeneous Poisson point process (PPP) with a unique receiver. The channel 
models span many blocks of time, and so the throughput variations over time can be characterized 
with our framework but not with prior frameworks. Thus, ETC may be better able to accurately 
suggest how to effectively use transmissions over time and space, such as multi-antenna transmission 
and opportunistic scheduling, and how much such techniques will improve area spectral efficiency 
in a long-term sense. 

A. Motivation and Related Work 

In the literature on wireless network throughput (see [HJ-[|3 and the references therein), a unified 
and time-invariant channel model is typically adopted over the entire network, but channels in a large- 
scale wireless network usually are diverse across time and/or space. Using a channel model without 
temporal correlation does not capture how channel states evolve over time and thus the impact on 
network throughput from the temporal (and spatial) variations of channels is ignored, and techniques 
which exploit the variations and correlations cannot be properly quantified. For example, we observe 
that transmission techniques that increase the snapshot network throughput may not increase ETC 
or may even degrade it. In particular, we will show that channel-aware opportunistic transmission 
(CAOT), i.e. transmitting when channels are in good states, does not always improve ETC, which is 
perhaps surprising. 

ETC requires the use of different channel models that include temporal discrepancies. We propose 
a finite-state Markov chain (FSMC) to model the fading channels, in particular a m-state Markov 
chain that is irreducible and positive recurrent. Each channel undergoes path loss and fading and has 
an ergodic property in that its fading state has an invariant (steady-state) probability flU. This idea 
can be traced back to the early work of Gilbert and Elliott [0 [TTOl that used a two-state Markov 
chain to represent good and bad channel conditions which was extended to a finite state case in [fTTTl . 

B. Contributions 

The first contribution in this paper is the model for ETC itself. We then calculate the ETC, which 
requires the outage probability for each fading channel state to be found, for which we find tight 
closed-form bounds based on a proposed 5-level interfering coverage area around a receiver. Any 
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single transmitter in the 5-level region of its unintended receiver will cause an outage if its interference 
power is enlarged by a factor 5. We show that appropriately choosing 5 admits bounds on outage 
probability that are much tighter than those found in previous work. 

Bounds on ETC and their corresponding scaling laws for some special cases are then found. They 
reveal several interesting implications, e.g. the scaling of ETC for a sparse or dense network is 



c E = e($ T s e ), (i) 



where C E denotes ETC, s f 
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in which §k is the function of the A;th state of 
a Markov fading channel with m states, $ = [cf>i 02 * • • 4>m] J £ [0, l] m is the invariant probability 
vectoiQ of the Markov fading channel model and a > 2 is the path loss exponent. From ©, we 
notice that a single deep fading state is not necessary to have a significant negative effect on ETC if 
its invariant probability is very small. This point is not revealed in prior work that neglects temporal 
variations. Also, we observe that ETC has a geometric interpretation because it can be viewed as an 
inner product of two vectors. Thus, ETC is maximized when the directions of and s e coincide. In 
addition, we show that channel-aware opportunisitc transmission (CAOT), which is a scheme to allow 
transmitters to transmit only when their channels are in good states, does not necessarily provide an 
ETC gain. Although CAOT is able to increase the transmission capacity contributed by good channel 
states, it loses the transmission capacity contributed by bad channel states. So CAOT cannot benefit 
ETC if the improvement is no larger than the loss from bad channel states. 

Three interference management methods - avoidance, suppression and cancellation - are applied to 
the network to reduce outage probability. Bounds on the outage probability and ETC with interference 
management are found which provide geometric insight into the efficiency of each technique in dif- 
ferent scenarios. For example, we show that interference cancellation is not effective for significantly 
increasing ETC in a spare or dense network. Also, we show interference management can control the 
direction and magnitude of vector s e . Finally, we show that CAOT should not be used if interference 
management can significantly lower interference. 

II. Network Model and Definitions 

A. The Network Model 

The wireless network considered in this paper is of infinite size and all nodes in the network 
are independently and randomly scattered. Thus, we employ a marked homogeneous Poisson point 



'The physical meaning of is the fraction of time that channel state sojourns as time goes to infinity. 
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process (PPP) II on the plane M 2 to represent the locations of all transmitting nodes in the network, 
which can be written at time r as^l 

n 4 {(X, G M 2 ,^(r) G R+,ei{r)) : e,(r) = 1} , (2) 

where X denotes node X as well as its location, H^t) is the fading channel gain between node 
Xi and its receiver and e^r) G {0, 1} represents the transmitting index of node X^. ej(r) = 1 
means node X{ is transmitting; otherwise, it is idle. The intensity (density) of II is A for all r G N. 

Each transmitter has a unique receiver and the distance between a TX-RX pair is a constant d > 1. 
All of the transmitters use the same transmit power and the channel model between each TX-RX 
pair is subject to path loss and fading. So the channel gain for TX-RX pair % can be written as 

H i (r)e(\X l -Y i \) = H l (r)£(d) J 

where Hi{r) is the fading channel gain, |Xj — Y$| denotes the Euclidean distance between nodes Xi 
and Yi and £(| • |) is the path loss function. In order to avoid the singularity where \X 0, we will 
use 

£(\X\) = G [1, oo)), X G M 2 , (3) 



where a > 2 is the path loss exponent and t(x G X) denotes the indicator function: t(x G X) = 1 
if x £ X and 0, otherwise. 

Specifically, we use an m-state FSMC model to characterize the fading effect of all channels in the 
network. The FSMC is irreducible and positive recurrent, and its m states are ordered. The FSMC 
model with transition matrix P is denoted by <S(P) G IR+ and S is an order set of the m states, i.e. 
for any two states Sj, Sj G S we have Sj < Sj where i < j and i, j G M. = {1, 2, • • • , m}. Since 5 
is irreducible and positive recurrent, the fading channel gain H(r) for all TX-RX pairs must satisfy 
the following conditions (H: 



<p k ^ lim }Vl(ff(r) G s fc ) and = V> fc = 1, (4) 

T = fc=l 

where $ = [0i,0 2 , • ' ' 5 0m] T i s m e invariant probability vector of P and u = [1 1 • • • 1] T is an 
m-tuple vector. Namely, at any time r, H{j) must belong to one of the states in S and <pk represents 
the probability that H{r) visits state Sk in a long-term sense. We can show this m-state Markov 
channel model has a temporal ergodic property as stated in the following lemma. 



2 Here it is better to use n(r) instead of II. However, to simplify notation, we will use II to stand for II(r) throughout this paper 
if ignoring time indices does not induce any ambiguity. This custom is applied to other set and variable notations. 
3 In a planar network, we require a > 2 to have bounded interference, i.e. It < oo almost surely if a > 2 1121 . 1131 . 
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Lemma 1 (Temporal Ergodicity of an instate Markov Chain): Suppose <S(P) is an irreducible 
and positive recurrent Markov chain with m states, and its transition matrix is P. Let h : S — > [0, 1] 
be a state measurable function of S and Z[t) is a random variable taking values in S. Thus, we 
have 



P 



1 ±j — _L 7ft 

lim - £ /»(Z(t) € 5) = X> /»(* 

L— >oo L z — * z — * 

r=0 i=l 



(5) 



where {</> fc , fc = 1, • • - , m} are the invariant (steady state) distribution of P and s k is the kth state 
of S. 

Proof: See Appendix |Al ■ 
The definition of ergodic transmission capacity in the following subsection is built based on the 

result in Lemma [T] In addition, the following lemma shows that the fading channel model of a FSMC 

also has a spatial ergodic property. 

Lemma 2 (Spatial Ergodicity): Consider a marked homogeneous PPP II with an independent 

mark H(r) E S, and let g : II — > M + be a measurable function on II. For any bounded subset 

A n C M 2 and fi(A n ) — > oo as n — > oo, we have 

E[<?(n fc )]= lim 7~a V f g((X,Hes k ))n(dX) = <f> k E[g(Il)}, a.s., (6) 

n,r->oo fl(An) J An 

where Il fe = {(Xj, Hi) e II : Hi e s^} is the PPP with channel state s k - 
Proof: Since II is homogeneous, we know 

Efo(II)]= lim -i- / E„[<7((X,i0)]/i(dX). 

Since {-f^i} are independent, 11^ is just the thinning homogeneous PPP of II and thus we have 

E[g(U k )} = lim — ^ / es fe )HdX) 

= lim — ±— / E[^((X,if))]E[l(iJG Sfc )]/i(dX) ( = ) fc E^(n)], 

where (*) follows from the temporal ergodicity result in Lemma Q] ■ 
Lemma |2l indicates that the spatial average of g(Hk) is equal to 0fcE[g(II)], So we know the intensity 
of n fc is Afc = 0fcA provided that g(-) is an intensity measure. 

The interference channel gain from transmitter Xj to its non-intended receiver Yi is denoted by 
Hji(r) i(\Xj — Yi\) where Hji(r) 6 S. The aggregate interference normalized by the transmit power 
at receiver Yi can thus be expressed as 

J f (r)= ]T ^(rKd^-y,!), (7) 
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where Ii is also called a spatial shot noise process |[T2l . lfl4l - lfT6l since it captures the cumulative 
effect at location Yi of a set of random shocks appearing at random locations Xj, and Hji£(\Xj — Yi\) 
can be viewed as the impulse function that gives the attenuation of the transmit power in space. In 
order to have a successful transmission for TX-RX pair i, the following signal-to-interference ratio 
(SIR) condition at receiver node must hold at time r: 

SIRi(r, A) = -^j-tK > /?, (8) 
d a Ii(r) 

where f3 is the SIR threshold for TX-RX pair i to successfully decode the received data. The network 
is assumed to be interference-limited. 

Note that according to Slivnyak's theorem [fTTI the statistics of Ii seen by any node in the network 
is the same if the nodes form a homogeneous PPR That means the average outage probability of 
each receiver node may be found by evaluating the SIR seen by a receiver located at the origin. 
Intuitively, the distribution of the point process is unaffected by the addition of a receiver at the 
origin, and this receiver is called a typical receiver. The performance measured at the origin is often 
referred to the Palm measure, and in keeping with simplified notation we will denote the probability 
and expectation of functionals of evaluated at the origin o by P and E, respectively. Also, Table [I] 
summaries the main mathematical notation used in this paper. 

B. Definitions 

Consider the typical TX-RX pair and its steady state outage probability is 

lim P[SIR(r, A) < /3] E {q k (X),k G M}, (9) 

T— >00 

where gfc(A) = lim T _ s . 00 P[SIR(r, A) < (5\H[r) = s&] is the outage probability for channel state Sk as 
time r goes to infinity. Now we are ready to use © to define ergodic transmission capacity in this 
paper. 

Definition 1 (Ergodic Transmission Capacity): Suppose transmitting nodes in a wireless ad hoc 
network form a homogeneous PPP of intensity A. For a given e E (0, 1), the ergodic transmission 
capacity (ETC) of a wireless ad hoc network is defined by 

C E ^6A E (l-e), (10) 

where b is the supportable transmission rate, e is the upper bound on the outage probability of each 
channel state, and A e = sup{A > : Ylk=i ^^(A) < e} is called maximum contention intensity 
achieved under the outage probability constraint e. 
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The definition of C E in (flOl) originates from the following definition: 



Ce = ^$>[SIR(t,A £ )>/3]. 



(11) 



r=0 



Since all channels are an irreducible and positive current Markov chain, according to Lemma[T|they all 
have temporal ergodicity. Thus, the definition in (fTTT ) is equivalent to C E = b\ E Y^k=i — <?fc(A E )]. 
This is the reason why we directly use the invariant probability of a Markov chain to define ETC 
instead of using (fTT|) . For ease of analysis, we need to quantitatively define the sparseness and 
denseness of a network with Poisson-distributed nodes. 

Definition 2 (Spatial Sparseness and Denseness of a Poisson-Distributed Network): Suppose the 
transmission coverage of a transmitter is the circular area of radius d. A network whose transmitting 
nodes form a homogeneous PPP of intensity A is called "dense" ("sparse") if the average number of 
transmitting nodes in the coverage is sufficiently large (small), i.e. nd 2 X 3> 1 (nd 2 <C 1). 



In this section, we study the general results of ETC. First, we have to calculate the outage 
probability for each channel state; however, only the bounds on the outage probability and ETC 
can be characterized due to the complicated distribution of the interference. According to the found 
bounds, the scaling behaviors of ETC are characterized and they reveal several observations. 

A. Bounds on the Outage Probability 



Since a closed-form expression of the outage probability defined in © is difficult to findj, we 
resort to bounds. The idea of approaching the lower bound for a receiver with channel state s k is to 
use a S-level interfering coverage for the typical receiver Y with fading state s k , and it is defined 
as follow 



which means any single interferer within it can cause outage at receiver Y with a SIR threshold 5(3. 
If Z| is not empty, it could contain dominant interferers and non-dominant interferers. In addition, 
U 5 k = (II nil) \ X is called 5-level interfering point process. 

4 If all channels are instead Rayleigh fading, the closed-form of the outage probability can be found by the Laplace transform of the 
aggregate interference contributed by Poisson-distributed transmitters [12] [18|. However, such closed-form outage probability cannot 
be obtained for the case of channels without fading |2] 1 19] or with a single state at any time. 

5 If 8 = 1, then Xf. is called the dominant interferer coverage in which a single interferer causes outage at receiver Yq. 



III. General Results on Ergodic Transmission Capacity 





(12) 



s 



The lower bound on the outage probability gfc(A) can be acquired by considering the outage events 
caused by U s k . The upper bound can be approached by finding the probability of the union outage 
events separately caused by the interferers in H s k and II \ nj*. These two bounds found are tighter 
than those in the previous works [|2) |fl9l , as the following theorem shows. 

Theorem 1: The outage probability ^(A) in © can be bounded as 



. 2 



l- e -H»> h a -*) < Qk ( X } < 1 _ (i _ A fc (A))+ e - A( ^ Q ~ 7r) , (13) 
where = max{x, 0} and A fc (-) is defined as 
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Afc(A) = ^ , (14) 

(s£ /d a 5l3-\ v ) 2 

and v, r], a 2 are respectively given by 

k=l K ' ^ k=l 

Proof: See Appendix |Bj ■ 
The physical meanings of v, r\ and cr 2 are the mean area of Z| with s\. = 1, the mean, and the 
variance of the interference contributed by the interferers of II \ 1T| for A = = 1. When the 
channel state Sk is high, the outage probability is reduced because SIR is large or equivalently the 
target SIR (3 is reduced. Nevertheless, it also can be explained from a geometric point of view. In 
([8]), we can let fading gain be incorporated into the path loss model of all interference channels, 
and according to the conservation property of a homogeneous PPP [17J, the intensity of the original 

2 _2 

PPP is changed from A to A/s£. This is why A in the bounds is scaled by s k a and thus interference 

2_ 

generated by the PPP with intensity A/s£ is small when Sk is large. So tightness of the bounds in 
(fT3l can also be observed. 

_ 2. 

If (•)+ in <fT3l is non-zero, the gap between the upper and lower bounds is A fe (A)e~ A( ^ s fc " which 
is a function of 8, A and s^. Since e~ u * , a and r] are all monotonically decreasing functions of 5, 

_ 2. _ 2. 

and the denominator of Afc(A)e~ A( ^fc " ~^ is convex for 5, it is easy to realize that A fc (A)e _A ^ s fc " ~^ 
is smaller than that without 5 if 5 is chosen appropriately. Fig. Q] shows the simulation results for 
channel fading modeled by a 2-state Markov chain. As expected, the two gaps decrease along with 
5 so that using 5 > 1 can make the bounds (much) tighter. In addition, the gap for a good channel 
state is much larger than that for a bad channel state. Hence, we should choose a sufficiently large 
5 in order to have tight bounds when the Markov chain has very good channel states. 

The result in (fT3l) will become slightly different if a transmitter uses a channel-aware opportunistic 
transmission (CAOT) policy. Recall that the states of a Markov fading channel are ordered so that a 
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better state has a higher subscript index. Suppose we call a channel state "good" in each FSMC if 
its subscript index is greater than or equal to g and s g > 1, which means channel gain H is good if 
H G S g = {s g , ■ ■ ■ , s m }. Therefore, the PPP with good channel states can be expressed as 

n g = {{X h Hi{r)) G II : Hi(r) G S g }. (15) 

According to Lemmas Q] and |2l its intensity is A g = XP[H G «S g ] = A Ylk=g $k as time goes to 
infinity. Therefore, the bounds on the outage probability with CAOT can be obtained from (TT~3T > by 
replacing A with A g , which yields 

2 2 

1 _ e -Ag("k S -T) < gfc(Ag) <!_(!_ A,(A g )) + e - x ^ s ^-*\ (16) 

Note that the bounds decrease in (IT3T ) when A decreases. Thus, the bonds in (fT6l decrease compared 
with (fT3l) . So CAOT improves the bounds on the outage probability of each channel state because 
nodes with bad channel states refrain from transmitting. In Section ITH-Cl however, we will point out 
that it may not always improve ETC. 



B. Ergodic Transmission Capacity 

By using the bounds on the outage probability of each channel state in Theorem [Q the ETC in 
(flOl) has bounds as shown in the following theorem. 

Theorem 2: Suppose the outage probability (j%(A) in © is upper bounded by e G (0, 1). Using the 
inequality in (fT3l) . bounds on the maximum contention intensity A E which maximize Y^k=i ^^(Ae) 
under the constraint of e can be given by 



J2 s lfaK< A E <-ln(l-e)]T 



(17) 



k=l 



where A| is given by in the following. 

At = : 



{A>0:iln 


[(1-A*(A))+- 


< 




1 - e 





k=i v - s£ir 



< us 



7T 



(18) 



where A fc (A) is defined in (fl4l) . 

Proof: See Appendix O ■ 
Scaling Laws of Ergodic Transmission Capacity. If we consider e — > 0, the upper bound in 
(fTTT) will reduce to eYILi Sk4>k/(v - s^tt) + 6(e 2 ) since -ln(l - e) -)■ e. So we know A^ : < 1 
and thus the solution in £[8]) is A^ = s% e/[(a 2 d 2a p 2 5 2 sl + " + v) - s^tc] + 6(e 2 ). We can conclude 
that the lower bounds in CLU) is e J2k=i <j>kS% /l(a 2 d 2a (3 2 5 2 sl + * + u) - s^tt] + 0(e 2 ) as e 0. So 
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the bounds in (fTTT) are asymptotically tight when the network is sparse because A = 0(e) as e — > 0. 
The scaling behavior of A E in this case turns out to be 

- Xt = e U±-*4-), d9) 

\ U k= i 1- TTS^/uJ 

where the argument in ©(•) only keeps the key parameters of interest. The above result is also the 
scaling law for a dense network because v\ <J — ln(l — e) as ird 2 \ is sufficiently large. Note that 
<n~9b is only valid for small A but large nd 2 \ since for small A we have to keep e small and nd 2 
sufficiently large in order to make nd 2 X large and the outage probability lower than e. In addition, 
if v — > oo (i.e. d and/or 5 — > oo), then (fT9l) further reduces to 

Ae = © ( — T I , (20) 

which means ETC is not affected by the fading channel states if the received signal is very weak. 
This makes sense in that channels can be equivalently viewed in a bad state all the time when the 
received signal is very weak due to long transmission distance. 

ETC with Channel-Aware Opportunistic Transmission (CAOT). The result in Theorem |2] is 
obtained without any transmission scheduling. Suppose now the channel state information (CSI) is 
available at each transmitter. Then transmitters can use their CSI to do CAOT and thus we have the 
following corollary. 

Corollary 1: If all transmitters transmit only when their channel fading gains are in S g , the bounds 
on the ergodic transmission capacity are 

^£4«<A E <^Mi^) E ^, (21) 

^9 fc=g ^9 k=qV-S%ll 

where <p Q = Y,k= q fa- 
Proof: For CAOT, the bounds on the outage probability for channel state Sk is shown in (fT6l) . 
Hence, the bounds in (I2TT) can be acquired by first taking off the terms with an index i lower than g 
in (fTTT) and replacing A E in the bounds with y? g A E . Since u, r] and a 2 do not depend on A g , we can 
replace A g in (fT6l) by A. Then following the same steps in the proof of Theorem |2] to derive (fT8l) . 
the lower bound in (1211 is completely achieved. ■ 

C. Observations and Discussion 

In the previous subsection, we have obtained bounds on ETC and discussed the scaling laws 
of ETC for a sparse and dense network. From the bounds and scaling laws, we have made three 
interesting observations. 
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ETC implicitly possesses a geometric interpretation. The scaling law of ETC in (fl~9l) can be 

expressed in a general form by vectors $ and s £ as 

C E = 6 ($ T s e ) , (22) 
where $ = [0i, </> 2 , ■ • • , m ] T and s is defined as 



A 

s f = e 



12. 2. 

q at Q a Q a 

oi Or) Om 



2_1 2_ 1 1 2_ 

V — TTS" V — 7TS% V — TTSm. 



(23) 



In other words, C E is scaled by the inner product of vectors $ and s e . So the result in (I22I) can 
be interpreted from a geometric point of view. Suppose the Markov channel model has two fading 
states (i.e. m = 2, this is so called Gilbert-Elliott channel model [0 [fTPl ). Then C E in (|22l) can be 
schematically presented in Fig. |2] Note that s e must be above on the 45° line because s 2 is larger than 
s 1 . The inner product of s e and $ can be written as $ T s £ = |$||s e | cos^ and 9 is the angle between 
vectors s e and $. So we will have a larger ETC if $ has the same direction as s e . The optimal $* 
that maximizes ETC can be given by = Therefore, if all Markov fading channels have the 
optimal distribution $*, then C E in (|22l) becomes 

and thus it is completely characterized by all channel fading states. 

Dominant channel states may not dominate ETC. Dominant states in a Markov chain means that 
their invariant probabilities are much larger than other states' invariant probabilities. In other words, 
if a channel has dominant states then it is in these states most of the time. Dominant channel states 
may not contribute too much ETC since their state magnitudes could be very small (very bad states). 
Thus, dominant channel states which really dominates ETC only when they have a large magnitude. 
This point can also be visually explained by Fig. [2l Suppose s 2 dominates and it is much larger than 
s 1 . In this case, $ and s e will move up and be close to the vertical axis. The projection of $ on s e 
will largely increase and it is mostly contributed by the vertical component. Thus, whether a channel 
state dominates ETC or not depends on the product of its magnitude and invariant probability. 

CAOT may not benefit ETC. If we compare the results in Theorem |2] and Corollary [Q we can 
find CAOT indeed increases the bounds with good channel states. Nevertheless, it may not always 
improve ETC since it loses the throughput contributed by bad channel states. To show this, let the 
upper bound in (fTTb be greater than the upper bound in (EH . This leads to the following inequality: 

i_i)f:^k T <g^k T . (25) 

^9 / k=gV — TTSj; k=l V ~ 1TS% 
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The LHS in the above expression is the ETC improved by CAOT and the RHS is the ETC loss 
because of bad channel states. If this inequality is valid, apparently CAOT may not improve ETC 
because the ETC increase in the good states could not compensate the ETC loss in the bad states. 
Hence, following from (|25l) . the better policy of using CAOT for a transmitter is when the following 
condition holds 



2 2 



^ ^ $^ • (26) 

However, the above condition is not implementable by transmitters in a wireless ad hoc network 
without knowing $ in advance. So CAOT is not always an effective means to enhance ETC in 
a real-time situation. A simulation example for ETC with and without CAOT is shown in Fig. [3] 
Channel fading is modeled by a 2-state Markov chain and the simulation condition is set to let the 
bad channel state be dominant. Obviously, we can see that ETC with CAOT is worse than ETC 
without CAOT. 

IV. Ergodic Transmission Capacity with Interference Management 

In Section ITlI-Cl we observed that refraining from transmitting when channels are in bad states does 
not necessarily increase ETC. This is because the CAOT scheme increases the transmission capacity 
for good channel states but further lowers the transmission capacity for users in bad channel states. 
That is, the throughput increase does not in general compensate for the throughput loss, particularly 
when bad channel states are dominant (i.e. channels are bad most of the time). Thus, the key to 
increasing ETC is to boost every entry of vector s e and not to sacrifice transmission opportunities 
of users with bad channel states. Two possible approaches to attaining this goal are through power 
control and interference management. In general, power control makes the analysis of the outage 
probability more intractable due to the complex structure of the interference. In addition, it is not a 
very effective means to increase SIR in a interference-limited ad hoc work lfl9l fl20ll . 

A. Interference Management - A Stochastic Geometry Perspective 

Interference management can be classified into three categories: interference avoidance, suppres- 
sion and cancellation. Avoiding interference in a wireless ad hoc network is typically accomplished 
by using space, time or frequency orthogonality to eliminate the co-reception of strong interferers. 
Frequency-hopping and CSMA are prominent examples of avoiding interference in an ad hoc network. 
Interference suppression deploys signal processing at the transmitter and/or receiver to (linearly) 
suppress interference without actually cancelling it. Direct-Sequence CDMA (DS-CDMA) is a typical 



13 



example of this category. In addition, receivers can try to cancel strong interference from their nearby 
unintended transmitters (e.g. successive interference cancellation (SIC) rt2TTl - [|23l ). Although avoiding 
interference and suppressing interference are two different methods of reducing interference, we know 
there exists a duality property between them in a wireless network with Poisson-distributed nodes. 
According to the conservation property of a homogeneous PPP [17 u, these two methods both reduce 



the original intensity of transmitters so that their effect can be demonstrated via another homogeneous 
PPP with a new intensity. However, interference avoidance has a better efficiency in reducing the 
intensity of interferers than interference suppression [0 [|20ll . 

The effect of interference cancellation can also be grasped from a geometric perspective. To explain 
this, suppose now any receiver Y { in the network is able to cancel some interference from its nearby 
interferers after the interference is avoided and/or suppressed. As time goes to infinity, the interference 
cancellation coverage of receiver Y with channel state s k is defined in the following. 

CI = \xeR 2 : > -J— \ , (27) 

where 7 is the interference of receiver Y and ^ k G (0, 1) is called interference reduction factor 
for channel state sJzl Coverage C£ means that any received interference within this region can be 
decoded by receiver Y with channel state s k , and all transmitters in C£ have a larger received power 
than transmitter X if C£ D (I1\X ) is not empty and f3 > 1. Also, C£ is a random compact set so that 
its mean Lebesgue measure /i(C£) is finite. If each receiver can perfectly cancel all interferers in its 
cancellation coverage after suppressing and/or avoiding some interference, then its SIR for channel 
state s k is 

SIR fe = S ^t , (28) 

£ Xjenr ^-(T)|X,K 

where il£ c = IT \ (C£ fl n U X ) is the noncancelable part of II with intensity X'jg . So equation (|28l) 
essentially suggests that interference management can be equivalently reflected by constructing a new 
PPP with a reduced intensity. The above SIR^ expression will be used in the following subsection 
to find the bounds on the outage probability. Those bounds are used to characterize the bounds on 
ETC with interference management. 

6 The conservation property of a homogeneous PPP with intensity A is that the intensity will change to A/a if all locations of the 
nodes in the PPP are scaled by a constant yfa. 

'Reduction factor 7^ can account for the joint effect of interference avoidance and suppression. For example, if there are M > 1 
available channels for DS-CDMA with spreading gain G > 1, then 7^ = 1/GM" 5 ". 
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B. Bounds on Outage Probability with Interference Management 

Bounds on the outage probability with interference management are shown in the following 
theorem. 

Theorem 3: If each receiver is able to avoid and/or suppress interference, and cancel interferers 
in its cancellation coverage, then bounds on the outage probability for channel state s k are given by 

1 - e-W"* "-*) < q k (X) < 1 - (l - A k (lh)) + e- x l^ % -*\ (29) 

where A£ = 7^ A ^1 — ^\ is the average intensity of the transmitters in X| with interference 
management and v\ is the mean Lebesgue measure of C%. 

Proof: See Appendix [Ql ■ 

There are a couple of observations that can be drawn from Theorem |3] First, the upper and 
lower bounds in (|29l are smaller than those in (TT~3T ). Canceling interference can be viewed as 
constructing a new PPP with reduced intensity 7^ A(l — Vj:/v) + in Z|, and it is more efficient 
to reduce interference than suppressing interference since it completely eliminates transmitters with 
strong interference power and thus the 1 — u^/u term does not have an exponent of -. Thus, we 
can infer that imperfect interference cancellation (i.e. interference suppression) is not as efficient as 
perfect interference cancellation (i.e. intensity reduction) and interference avoidance since it merely 
decreases transmitters' interference and does not directly reduce transmitter intensity. 

Second, interference cancellation is not equally useful for all networks. For example, in a dense 
network, canceling the strong interferences from the nearby transmitters can significantly reduce 
outage probability such that network throughput is substantially increased. This point can be easily 
verified by letting \v be sufficiently large. In this case, qk{X) is close to unity if no interference 
is canceled. On the other hand, for sparse networks, interference cancellation may merely have 
a marginal reduction in outage probability. For sufficiently small Xv and v > vZ, (l29l) can be 

2_ 

simplified as qk(X) = ( a X(u — v^) + 0((Xu) 2 ). So when receivers cancel more interferers, its 
outage probability is reduced by amount of 0(Xu k ) which is really a small and trivial improvement. 
From this observation, we see that canceling strong interferers for each channel state may not be an 
effective means to increase transmission capacity for small e since the maximum contention intensity 
of transmission capacity for each channel state is already a small value in this case. 

C. ETC with Interference Management 

According to the bounds on the outage probability in Theorem [3l ETC with interference manage- 
ment is bounded as shown in the following theorem. 
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Corollary 2: Suppose interference management is used in the network and each receiver is able 
to perfectly cancel all interferers in the interference cancellation coverage of each channel state. Let 
the outage probability be upper bounded by e G (0, 1). If C£ c Z|, the maximum contention intensity 
for channel state has the bounds given by 

£ ( - ) " « < A E < - ln(l -e)J2 (-) " , (30) 

where v > v% and X% is given by 

XI = inf J A > : iln 



(1-A*(A))- 



1 - e 



<\ l -i)[v- s k*)\- (3D 



However, if Xl C Ci , then (l30l) becomes 



A E > ( - ) " ^> (32) 



where A| = inf |a > : (l/d a 6(3 - Xrjsl l lk < 4 "7^t A }- 

Proof: See Appendix El ■ 
By comparing (|3Q|) with (TTTI) . we can perceive that the effect of interference cancellation on 
ergodic transmission capacity can be interpreted to shrink v by (1 — z/£/V)-fold. This is equivalent 
to saying channel gain s k increases (1 — v%jv) ~t -fold. Suppose e is small and thus from (|30l and 
(I3TI) we know A E ~ ^Y^k=i(~)" ■ ^° interference cancellation in a sparse network 

Tfe (i/-irs£)(l-j/£/i/) 

can make the transmission capacity for channel state increase (1 — z/^/z/) _1 -fold. Although 
avoiding and suppressing interference can linearly augment ETC in a sparse network, interference 
cancellation could contribute much more ETC than them if v k is very close to v. Fig. |4] presents 
a simulation example showing how interference management improves ETC. We first notice that 
interference management does not provide too much ETC gain when e is extremely small. On 
the other hand, if the network is very dense, the efficacy of interference cancellation is seriously 
weakened because interference is large. So the solid-circle curve of ETC looks like a concave function 
of e. Therefore, interference management in an extremely sparse or dense network can merely have 
marginal improvement on outage probability. 

How should interference management be used for each channel state to maximize ETC? In Section 
IIII-Cl we have pointed out that ETC has a geometric interpretation since its bounds can be viewed 
as the inner product of two vectors: vectors $ and s e should roughly align. Since vector $ is a 
channel characteristic, it cannot be manipulated to the desired direction. Therefore, the only option 
is to design vector s e such that it is enlarged and rotated to the direction of $ as closely as possible. 
This can be attained by interference management. To illustrate the idea of how to change s e , the 
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r 1 2.1 T a f 2\~ 

right part of Fig. |2]is redrawn in Fig. |5] Let s e = s" s£ where §k = ( 

2 2 

s * = ( s i)" ( s 2)" represent the optimal vector that s e can achieve by interference management. 

2_ 

Note that s* k = — (1 — is£/v)~%, k = 1,2. Therefore, after interference is reduced, vertices c (sf , 0) 

— 2 2 

and a (0, s% ) can be maximally pushed out to vertices h ((s*)«, 0) and f (0, (s^)«). Vertex i is the 
point where s* is projected on $ and Vertex d is the point where s e is projected on $. The distance 
from m to e represents the increase of ETC due to interference management. 

Since s* is the best vector s e can achieve, how can we make vector s e move to vector s*? Namely, 
how should we choose 7^ and v% for each channel state such that s t can approach s*? The policy 
is to reduce interference for each channel state as much as possible because we can formulate the 
following nonlinear programming problem to optimize all 7^: 



v — 7is% ) and 



maX 7fc 



2 — 

f'fJaLA' ! r (33) 

tlX*-™;) 7* [(i - " l °/")] T 

subject to 7fc>7min fc , for all k G M, (34) 

where 7 m i nfc is the lower bound of 7^ and it can be determined by the system resources or limitations 
such as number of available channels and the maximum spreading gain, etc. Note that u% is a mono- 

CM 

tonically decreasing and nonlinear function of 7^ so that l/jk [(1 — v^/v)] 2 is also a monotonically 
decreasing function of 7^. Therefore, the optimal solution of 7^ must happen at 7^ = 7 m i nfc , which 
means interference should be avoided, suppressed and cancelled as much as possible in order to 
achieve s*. In addition, using interference management could make CAOT perform poorly because 
it may make most of channel states become "good" so that ip g fa 1. 

V. Conclusions 

In this paper, we presented a long-term look at the transmission capacity problem, which is 
completely different from the previous works on investigating network throughput at a particular 
time point. The motivation of this work is to understand how the temporal characteristic of a channel 
influences the network throughput with an outage probability constraint. Therefore, all channels are 
modeled by a m-state FSMC that has temporal and spatial ergodic properties. Bounds the on outage 
probability of each channel state and ETC for the case with and without interference management 
are all found and they show that ETC can be characterized by the inner product of vectors $ and s e . 
For a sparse or dense network, the scaling law of ETC that is derived from those bounds provides 
some guidelines on when to use CAOT and how to do interference management. 
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Appendix 

A. Proof of Lemma [7J 

Proof: We need to show that j;^([fi(2(r)) converges to Y1T= 1 ^ M s *) as L — )• oo almost 
surely. Suppose {Z(r), r > 0} has an invariant distribution {</>j, z = 1, 2, • • • , m} and define V^(L) = 



St=o is the number of visits to state s, before L. Since ft is a positive function, for any 



J C 5 we have 

L-l 



h(si 



L 

Vi(L) 



h( Si ) + ^ 



L 



L 

Vi(L) 



h(si 



L 



L 



ft(si) + 2 



where (*) follows from Pflirn^oo Vi(L)/L — fa] = 1 for all z = 1, 2, • • • , m. 

For a given £ > 0, choose with an appropriate size and consider L is sufficiently large so that 



^2faH(si)<- and 
Therefore, when L is sufficiently large it follows that 

L-l 



Vj{L) 
L 



< e, 



r=0 i=l 

which establishes the desired convergence. The proof is complete. ■ 
B. Proof of Theorem UJ 

First of all, we have to find the intensity \§ of ilf . According to [|24| . the Laplace functional of 



a homogeneous PPP IT for a nonnegative function w : 



is given by 



C n (w) 



E 



-f R2 w(x)n(dX) 



exp 



A(l 



-w(X) 



)n{dX) . 



(35) 



Since the Laplace functional completely characterizes the distribution of a point process, we can find 
the intensity of IT^. by calculating C u s(w). For a bounded Borel set A C M 2 , The Laplace functional 
of U s k with w(X) = w(X)l n s(X) can be written as follows: 



C n s(w) 



-Xfj,(A) 



E 

8=0 

oo 



A' 



A 



no 



-w(Xj) 



F[Xj e U s k ] + F[Xj i n*]) f x(dX 1 ) ■ ..p(dXi 



Ell/ (e~ 9[Y) nY e nj] + 1 - P[K e nfl) A/i(dY) 



(a) 



cxp 



i=0 

-A / (1-e 



8=1 



Si e 



s k d~ 



l{\Y\)Sfi 



oo 



MdE) 
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where (a) follows from the property of spatial ergodicity. Letting A — > M 2 and according to (1351) . 
we know the intensity of n£ is 

s k d~ 



i=l 



J(x)8/3 



oo 



x G 



(36) 



So U s k is a non-homogeneous PPP since X 5 k depends on x. 



Since il|(Z|) is a Poisson random variable, its mean can be found as follows: 



E[U 5 k (l s k )\ = E 



Xjeu\x 



(b) 



E[t ui (X)UdX) 



X s (\X\)n(dX) = 2nX 



i=l 



Si e 



s k d~ a 
J(x)5/3 



oo 



xdx 



2ttA^. 



i=l 



xdx = ttXs^ d 2 {5(3)' ^fas? -sM= \(s k °v-ic). 



i=l 



where (b) follows from the Campbell theorem iflTl . Let S{Ii 5 k ) denote the outage event caused by 
any transmitters in IT^ and its probability is 

P[5(n*)] = 1 - exp (-A (vsf« -*))< ft(A), 

which is a lower bound of (&(A) because it ignores the interference contributed by the transmitters 
that are not in n|. 

Let S°{U s k ) be the complement event of £(U s k ) and £°(Il s k ) means the outage event caused by the 
transmitters of II \ U s k . So the upper bound of qk(X) is given by 

g*(A) < P[£(n£) u s c (u s k )} = 1 - e~ x ^ § -^ + p^nf)^-^ 1 -^, (37) 

where P[£ c (n^.)] = F[s k d~ a < [31%] and I k is the interference contributed by the transmitters of 
II \ 11^. Unfortunately, it is impossible to explicitly calculate P[£(Z£)] and thus we resort to find 
its upper bound by Chebyshev's inequality. Using Campbell's theorem, the mean and variance of 



interference I k can be calculated as follows: 



E[I, 



E 



Ml*il)ln\n{(*J 



E[H] / £(\X\)[\-\ 5 (\X\)}^dX) = \sl 



Var[4 c ] = E [(&] - E[I c k r = E[H 2 } / [£(\X\)} 2 [A - X S (\X\)] fi(dX) = AaJ" o\ 

where a 2 is bounded is due to bounded r\ since [^(x)] 2 < £{x). The upper bound of P[£ c (n^)] can 
be obtained by 



P[£ c (n£)] = F[d~ a < 5(311} < P 



d~ a 
\d~ a 



s k «6{3X V < 5(3\P k 



5fiXr]\ 



\dr 



5f3Xr)\ 



< 



/d a S(3 - At?) 5 
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where (c) follows from Chebyshev's inequality. Substituting the above result into (|37l) . the proof is 
complete. 



C. Proof of Theorem [2] 

Since Y^k=i 0fc<?fc(A) < e, the lower bound in (fTJt with A E can be rewritten as 

m m 

^0 fc exp (^-X E (us k a - 7r)j > 1 - ^(pkQkM > 1 -e, 

which gives us 



fc=i 



fc=i 



fc=i 



-A E (ys, 



k=l 



k=l 



where At = argmax A Ae A (" s * a n \ Hence, e A fc^ s fc Q ^ > 1 — e so that we have 



m(l ~ e) , = h(l z f . § 



(38) 



S h "l/ — 7T fc=l 



So the upper bound on A E is acquired. 

Similarly, we know 1 — Yuk=i 'PkQkQ^v) > 1 — e and (fT3l) with A E can give us another lower bound 
on 1 — Y^T=i 'Pklki^s)- Combining these two lower bounds, it yields the following result: 



1 - ^2<piqi(X E ) > max ^ 1 - e, ^ (f> k (l - A fc (A E )) + e~ 



■A E (^ fe °>-w) 



(39) 



k=l 



k=l 



Since (1 — Afc(A E )) + e XE ^ us k a "O is a monotonically decreasing function of A E , the lower bound 
on A E must happen when (1 — Afc(A E )) + e~ AE( ^ fe is equal to 1 — e, which means A| = inf{A : 

_ 2. _ 2. 

(l-A fc (A)) + e- A ( 1/s fe a -^ < 1-e}. We can explicitly write the relationship (l-A fc (A))+e- A ^ a ~^ < 
1 — e in the following. 



In 



:i-A fc (A)y 

fl-e) 



which implies 



A; mHA>0:-j-m 



1 - e 



< ( I/S fe a - 7T 



Thus, the lower bound on A E can be written as ^2™ =1 4>kS% X%. The proof is complete. 
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D. Proof of Theorem \3\ 

Since the interference generated by II is scaled by j k , it is equivalent to the interference generated 

_2 

by a homogeneous PPP with intensity 7^ A. The intensity of il£ c at location X can be shown to be 

H\X\~ a /3 



A£ c (|X|)= 7 fAP 



The average number of nodes in C k is f^°[j k ^ — A£ c (|X|)]/i(dX) = 7^ \v k . The lower bound on 
the outage probability can be characterized by the average number of noncancelable nodes within 

_2 

2|, that means we have to find \[i(C k HX|) which is the average number of nodes in II£ C HX|. 
Transmitters in Xf and C£ must satisfy the following in two equalities, respectively: 



5(3Hj 

Sk 



l/a 



d (for Xj el s k ) and < 



P(d a lk Io + s k ) 



l/a 



d (for X 3 eCl). 



So for any X, e If. fl C|, we must have 

\Xj\ < d I mm 



Sk 



1 + 13 
' 6(3*(l + d° lk Io/sk) 



l/a 



(40) 



If C£ C Xf a.s., then we must have 



5 > 



l + (3 



a.s. 



5 > 



1 + /3 



^{l + d^kh/skY ' ' " (3 2 ' 

2_ 

In this case, the average number of noncancelable 5-level interferers is \^^{jj,{lf) — fJ.(C k )) which 

2 

is equal to \^£v{l — v k l v )- ^ n me °ther hand, if If C C£ a.s., 5 < ^-0- and thus z/£ > v in this 
case and thus all interferes in if are cancelable. Combining these two cases, the average number of 
noncancelable 5-level interferers should be written as z/A£. Using A£ to replace A of the lower bound 
in (TT3T) . we can have the lower bound on the outage probability. 

The (-) + term of the upper bound is due to the outage caused by the transmitters out of if so 
that the intensity of II \ IT| is 7^°A because no interferers are canceled in II \ IT^. So we can just 
replace A in the (■) + term of (IT3T > by 7^° A and also replace A in the exponential term with X k . Then 
the upper bound is obtained. 



E. Proof of Corollary \2\ 

By considering the given condition Y^k=i 0fc9fe(^E) < e, the success probability for channel state 
Sk obtained from (|29l ) is bounded as follows. 



1 - ^ 0fc£ < 1 - ^ 0fc<?fc(A E ) < ' 



— 2/ol . m 



fc=i 



k=l 



k=l 
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If C£ c Zf, then using the above inequality and following the same steps of finding the upper and 
lower bounds in the proof of Theorem |2] we can show the results in (|3Q|) and (13TI) . If if C C£, then 
the upper bound on success probability is no longer available since all 5-level interferers are canceled 
and thus [y — z/£) + = 0. As shown in the proof of Theorem |2] there are two lower bounds on the 
success probability: one is 1 — e, the other is Ylk=i $ k ~ ^ k (^*7fc )) ■ Considering A fc (-) < 1, 
the lower bound on the success probability can be expressed as 



in I m ^ 



fc=l k k=l 

which yields the following inequality 



? x ftZd<*8p-Xn 

_1 / i_2. _2 \ 2 a 1 

A > : s£ 7 fc a ( - Ar^s^ a 7 fe Q J < A^- >, 
which renders us the lower bound Y^k=i( s k/lk) "0fcA|. This completes the proof. 
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TABLE I 

Summary of Main Mathematical Notation 



Symbol 


Definition 


n 


Homogeneous PPP of transmitters 


A 


Intensity (density) of JT 


C E 


Ergodic transmission capacity 


A E 


Maximum Contention Intensity 


e 


Upper bound of outage probability 


d 


Transmission distance of a TX-RX pair 


S G R+ 


m-state Markov chain for modeling fading 


Sk 


fcth state of Markov chain 5 


4>k 


Invariant (steady state) probability of channel state Sk 


H(t) 


Fading channel gain at time r, H(t) G 5 


a > 2 


Path loss exponent 




SIR threshold for successful decoding 


n 


5-level interfering coverage for channel state st 


5 > 1 


Parameter of defining if 


V 


Mean area of if for Sk = 1 




Interference cancellation coverage for channel state Sk 


Ik G (0, 1) 


Interference reduction factor for channel state Sk 




Path loss function 


9fc(0 


Outage probability for channel state Sk 




Lebesgue measure of set A 
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Fig. 1. The gap between the upper and lower bounds on gfc(A). Channel fading is modeled by a Markov chain with 2 states. The 
network parameters for simulation are: d = 5m, A = 0.01, a = 3, /3 = 2, si = 0.5, S2 = 2 and cj>\ — <j>2 = 0.5. 



2. 
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Fig. 2. The Gilbert-Elliott channel model and its corresponding geometric presentation of Ce, where {pij} are the state transition 
probabilities for the FSMC model and Sk = Sfc/(^ — tts£ ) "2 . 
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Fig. 3. Numerical results for ETC with and without CAOT. The network parameters for simulation are : d — 10m, f3 = 2, a = 3, 
8 = 1.5, si = 0.5, s 2 = 2, </>i = 0.8, <j> 2 = 0.2 and A = 0.01. 




Upper Bound of Outage Probability e 

Fig. 4. Numerical results for ETC with and without interference management. The network parameters for simulation are : d = 10m, 
P = 2, a = 3, 8 = 2, si = 0.5, s 2 = 2, 01 = 2 = 0.5, 71 = 72 = 0.6 and A = 0.02. 
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Fig. 5. The geometric representation of the ETC for an FSMC with 2 states, s* is the optimal vector that s 6 can achieve by 
interference management. By using interference management, vertices a and c can be optimally moved to f and h, respectively. The 
projection points of vertices b and g on $ are d and i, respectively. 



